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Abstract 

' 

In this paper, we present a new variation of dilated matrix inequalities (Mis) for Bounded 
' Real MI, invariant set MI and constraint MI, for both state and output feedback synthesis 

, problems. In these dilated Mis, system matrices are separated from Lyapunov matrices to 

' allow the use of different Lyapunov matrices in multi-objective and robust problems. To 

D , demonstrate the benefit of these new dilated Mis over conventional ones, they are used 

in solving controller synthesis problem for systems with bounded actuator in disturbance 
l/^ , attenuation. It is shown that for the resulting multi-objective saturation problem, the new 

form of dilated Mis achieves an upper bound for L2 gain that is less than or equal to the 
upper bound estimate achieved by conventional method. 

: 1 INTRODUCTION 

■ Most of the (Linear) Matrix Inequality ((L)MI) characterizations in control techniques such as 

Hoa, H2, use a quadratic Lyapunov function {V = x^Px |P > ) to develop their Mis (e.g., 
see ref [1] ) . The resulting Mis end up having entries with the products of the Lyapunov variables 
and system matrices. This causes some degree of conservatism in multi-objective and robust 
J> . problems by forcing common Lyapunov matrices for all objectives. For example, see Ref. I2]-11] 

which make use of common Lyapunov variable in their multi-objective problems. 



(N 

^ . 

, Recently, researchers have been using matrix dilation results to reduce this conservatism. 

The earlier, and the best known, of the results obtained by these new techniques are in discrete 



. time settings ([6], [7]). Although a lot of effort has also been made for the continuous-time case, 

it is still an open problem, mostly due to the fact that dilations to reduce conservatism destroy 
the convexity in some important cases. Some of the very nice and convex results in continuous 
time are achieved by Ebihara et al. in [8]-|10j. in case where the problem can be case as 

X ■ AX + XA^ + 5iX + 5^AXA^ + XA^AX < (1) 

for a suitable choice of X, A, etc. By assigning different matrices to A, Si, 62 and A, this 
general form covers some continuous-time control problems such as stability, H2 and D-stability. 
Furthermore, references [15j and [16] present a technique that, using a projection (or Elimination 
lemma) based approach, leads to set of convex search for several important problems; H2, 
stability, D-stability, etc. 

Unfortunately, neither of the two approaches above deal with synthesis inequalities faced 
in the L2-gain (i.e., bounded real problems) or several invariant set determination problems. 
Roughly, these problems result in a constant term in ([1]). There have been dilated Mis for such 
cases, as well (|11]-[T3]. among others). So far, there seems to be two weaknesses associated 
with these set of results. Often, the synthesis results are for the full state case, by exploiting a 
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structure that holds only in full state problem. Furthermore, they all seem to need an additional 
scalar variable which enters the MI in a non-convex fashion ([Il]-[l3]). While it seems that such 
a non-convexity is inevitable, it can be addressed with a line search since the culprit is a scalar 
variable. 

Here, we present a new dilated MI, that can be applied to the bounded real LMI, as well as to 
the matrix inequalities that are used in the invariant set for peak bounded disturbance ([!]) and 
the constraint LMI ([!]). These new Mis are obtained explicitly through a constructive methods, 
to avoid the ambiguities that can - at time - accompany results based on the Projection Lemma 
(see [8] for a discussion on this issue). We also rely on a scalar variable that renders the 
problem non-convex and use line searches to obtain the final result. Fortunately, we show that 
the proposed approach is rather easily extended to the output feedback problem, assuming the 
controller is of the same order as the plant. 

Finally, we study the effect of this new matrix dilation technique in reducing conservatism 
in controller design for linear systems with bounded actuators which avoid saturation. As 
mentioned in [2], this is often set as a multi-objective problem, and can suffer the conservatism 
forced by common Lyapunov matrix. We show that, these new Mis lead us to a problem whose 
performance is at least equal to the one with standard Mis. 

The system we study has the standard model 

Xp — Axp + Biw + B2U 

z = CiXp + Diiw + D12U (2) 
y = C2Xp + D21W 

with the closed loop 

( X = Acix + Bciw 

\ z = CclX + DclW 

where the details differ in the state feedback and output feedback cases. The Transfer function of 
this system is Tzw{s) = Cci{sl — Aci)~^Bci + Dd- We use He(^) as short notation for A + 
to save space. The rest of the notations throughout the paper follow standard practices. 



2 Dilated Mis for Some Practical Design Specifications 

In this section, we derive equivalent dilated Mis for standard Bounded real LMI, invariant set 
Ml and constraint LMI. In the new Mis, the system matrices and the Lyapunov variable are 
decoupled. 



2.1 A Dilated LMI for L2 Gain 

Lemma 1 (Bounded Real Lemma p]) Ad is stable with \\Tzw{s)\\oo < Icon if cind only if 
a{Dci) < jcon CLnd there exists Qi > such that 

AaiQi + QiAl B,i QiCl 

* -jaoni Dl I < (4) 



Theorem 1 The closed-loop system ^ is stable and its L2 gain is less than jnew if there exist 
a positive constant < ei < 1, and square matrices Xi > and Gi which satisfy: 



( ^1 Bel 

-k ★ 



-Xi\ 





+ He(Q^GiP) < 



(5) 



where 



V 



[I - 26i/]; Q = [{Al--I) Cl I] 



Proof 1 Proof is through showing Mis ^ and ([^ are equivalent. Suppose that MI ^ holds. 
Consider the explicit bases of nullspaces of V and Q 
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/ / 0\ 
0/0 
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V-A^ + i/ -CJJ 



By multiplying Mq and its transpose from the right and left sides respectively, and considering 
QMq = 0, inequality (5\) becomes 



AciXi + XiA^i Bel ^1^1 



-k 
-k 



cl 

~^newl 



< 



(6) 



Now if we define Xi = Qi and "fnew = 7con; it is clear that ^ holds. Next, multiplying ^ 
from right and left by M-p and its transpose respectively gives 



Bel 







< 



(7) 



This inequality implies that (0j can have a solution only for < ei < 1. 

On the other hand, suppose that holds with Qi > 0. Note that ^ can he rewritten as 
(0j by defining Qi = Xi and ^con = Inew Since Xi > 0, for any ei > we have 



7^^(4eiXl)-l7^ > 



(8) 



where IZ = [— 2eiXi(A^ — |) — 2eiXiCj]. Since the right hand side of the inequality is of 
order ei, it is possible to find a sufficiently small ei > which for this ei, and any ei < el, the 
following holds 

[left side of 0;] + [left side of fS)] < 0. 
Applying the Schur complement to this inequality leads to 



B, 



I 

-k -k 
■k -k 



Del 

~^newl 
■k 





'2eiCciXi 
-4eiXi 



< 



where H = He((^c« ~ + ^i- choosing Gi = Gj = Xi, this inequality can be written 

as i.e., satisfaction of leads to a specific choice for the matrices Gi and Xi that satisfy 
1^, with the same performance estimate. 



2.2 A Dilated MI for the Invariant set 

Lemma 2 ( Invariant Set MI [Ij) £ = {x\x'^Px < w"^^^} is a reachable set (invariant set) 
for the LTI system ^ exposed to a peak bounded disturbance w{t)'^w{t) < i/^^ax there exist 
a scalar aeon > and Q2 = P^^ > such that the following MI is feasible 

f AelQ2 + Q2A^i + acouQ Bel 




Theorem 2 Inequality (0j is feasible if and only if there exist a square matrix G2, a constant 
< €2 < 1 and X2 = P^^ > which satisfy the following condition for some anew > 0/ 

X2 Bel — -^2 \ 

* -a„e»/ + He(Q^G2P) < (10) 

* * / 

where V = [I - 2e2l] and Q = [{A^ - i/) /], with A = A^i + ^I. 

Proof 2 Proof is similar to the Theorem 1. In this case TZ in inequality ^ is 

n=[-2e-2X2{A^ -h) 0] 

2.3 A Dilated MI for Constraint LMI 

Constraint LMI ([Ij): Suppose that y = Kx and x € {x\x^ Px < w"^^^}. Then, ||y|p will be 
less than or equal to uf-^ if the following LMI holds for Q3 = P"^ 

Theorem 3 The matrix inequality is feasible if and only if the following inequality is 

feasible for some < €3 < 1 and X3 > 0; 

^3 -X3 

* I + He(Q^G3P) < (12) 

* 

where V = [I - 263/] and Q=[-I - I]. 

Proof 3 Proof is similar to the Theorem 1. In this case TZ in inequality ^B^ is 

7e= [26-3X3 2e-3X3i^^] 

3 A multi-objective problem solved using the new dilated Mis 

In this section, the new Mis are used to solve a multi-objective problem. The problem considered 
here is that of controller design for a system with bounded actuators exposed to a peak bound 
disturbance, w{t)'^w{t) < w^ax- Nevertheless, the idea can be applied to other multi-objective 
cases. The goal is to design a controller which makes this system internally stable and guarantees 
disturbance attenuation, while avoiding any violation of saturation limit (uum)- As mentioned in 
[2], this problem is a multi-objective problem, often yielding very conservative results. Typically, 
the linear or low gain controller for this problem is based on finding a controller (a state feedback 
or dynamic output feedback compensator) which gives the best L2 performance while making 
sure that saturation limits are not violated, by keeping the maximum of the control input below 
the limit for all the points in the reachable set of the closed loop system. We can formulate the 
proposed solution through the following two algorithms: 



Algorithm 1 (Conventional approach): For common Qi = Q2 = Q3 = Q > 0, minimize 
7co„ in LMI dH) subject to Mis (l9D and (fTTTl . 



• Algorithm 2 (New approach using new dilated Mis): For Xi > 0, X2 = X3 > and 
common Gi = G2 = G3, minimize jnew in © subject to Mis (fTO|) and (fT2]l . 

In Algorithm 2, as we show later in both state feedback and full order output feedback 
cases, common G is needed to turn the Mis into appropriate form to get a unique controller. 
However, there is no obligation to use a common Lyapunov matrices for the L2 and invariant set 
inequalities. The reason to take X2 = A3 in Algorithm 2 is that we are keeping the maximum 
of the controller, below its limit, for all the points in the reachable set of the closed-loop system. 
Here, the reachable set is x-^A^^x < w^ax- 

The following theorem states the advantage of the new approach expressed through Algo- 
rithm 2 over the conventional one obtained by Algorithm 1, by guaranteing better (or at least 
no worse) L2 gain estimate. 

Theorem 4 (Multi-objective) For the multi- objective saturation problem mentioned above, 
Algorithm 2 with a common auxiliary variable, G, but with non-common Lyapunov variables, 
always achieves an upper bound estimate for the L2 gain that is less than or equal to L2 gain 
performance estimate achieved by Algorithm 1. 



Proof 4 If Algorithm 1 is solved, then Theorem [7] implies that there is a positive e such that 
for any ei < e by taking Xi = G = G^ = Q and jnew = Icon, we can satisfy ^ with the same 
closed-loop system derived by solving Algorithm 1. Similarly, following the same argument. 
Theorems 2 and 3 guarantee that there are €2 and £3 such that X2 = G = G^ = Q and 
Oinew = Oicon Satisfy ilO\) and il2\). Therefore, all the Mis in Algorithm 2 are feasible with 
Inew = Icon if wc sct G = G^ = Xi = X2 = Q and use the same closed-loop system, as 
obtained by Algorithm 1. 

Therefore, any solution of Algorithm 1 can be achieved by Algorithm 2, for a small enough 
e 's, without exploiting the ability to use different Lyapunov matrices, which could only improve 
the results. 

Recall that if MI ([5]) holds for some ei, it would hold for any e < ei. Same argument is true 
for Mis (fTOj) and (fT2]) . Therefore, using the same e still allows the results to be at least as good 
as those from Algorithm 1. Therefore in Algorithm 2, to decrease the computational cost of 
line search for ei, 62, and £3, we can use the same e for all Mis. This leads to some degree of 
conservatism. The conservatism is due to the fact that the best result obtained by Algorithm 2 
is not necessarily preserved if we lower all or any of the e^. For best results, the e's are allowed 
to vary independently (see the numerical examples). 

In Section 13.11 Algorithms 1 and 2 are applied to state-feedback case and in Section 13. 3| 
these algorithms are used for full order dynamic output feedback compensator. The result of 
these two algorithms are compared through numerical examples. 

3.1 State Feedback Case 

In this section, state feedback controller synthesis is considered. Therefore, u = Kxp, and 
the closed-loop system in (3) can easily be obtained with x = Xp. For state feedback control. 
Algorithm 1 and 2 can be expressed in the convex MI set as in lemmas below. 

Lemma 3 ([2]) System ^ with w{t)'^w{t) < w^^^. and state feedback K is internally stable, 
never saturates and has a disturbance attenuation level 'ycon 

if there exist Q > 0, Y and a 

positive constant Ocon > such that 



-Ico'rJ 
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n 



^conQ 



Bi 



< 



-Q 



-I 



< 



where H = AQ + QA^ + B2Y + B2- The variables in this problem are Q, , Y, "jcnn, and 
Oicon where aeon is searched through a line search. The controller is given by K = YQ~ 



icon ! 
-1 



Lemma 4 System ^ with w{t)'^vu{t) < tf^ax '^'^'^ state feedback K is internally stable, never 
saturates and has a disturbance attenuation level ^new if there exist Xi > 0, X2 > 0, Y , square 
matrix G, constant anew > 0, and small positive scalars < 1 (i = 1,2,3) such that 



/Xi + n + n'^ Bi 



•k 
•k 



-k 



12 



-Xi + G^-2ei(n)\ 


-2eiiCG + D12Y) 
-2ei{G^ + G) 



< 



'X2 + n + n^ + ii^(G + G^) Bi 

•k ★ 



-X2 + G^ - 262(11 +^Gy 


-262(G^ + G) 



< 



/X2 -G-G^ 



\ 



-Y 



-k 



G'^-2e3G' 
2€^Y 
-2e3(G^ + G) 



< 



where 11 = AG + B2Y — hG. The variables in this problem are: Lyapunov matrices Xi and 



X2 as well as jnew, Y, Oinew o,nd ei [i = 1,2,3) where ar^ 
search. The controller is given by K = YG~^. 



and e 's are searched through line 



In this lemma, same slack variable G (Gi = G2 = G3 = G) is used to make the MI set 
convex and to get a unique solution for K. 



3.2 Numerical Example: State Feedback 

In this section we use the results of Section 13.11 in a numerical example. Consider the example 
from [13] 
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with k = 2 and / = 0.2. The system is exposed to a peak bounded disturbance with Wmax = 5. 
The controller limit is = 8. Using Algorithm 1, we get the following controller 

K = [ -1.7970 -0.7094 -2.2916 -2.1091 ] 

which gives the minimum 7*^^ = 1.3038. This problem is also solved through Algorithm 2. 
To decrease the computational cost, we used a same e in all three dilated Mis. Based on the 
simulations that we have done so far, the variation of ^new with e is bowl shaped. Therefore, 
instead of a full line search on e, we used the Golden-ratio method to identify the minimum 
performance, ")newi ^'^d the corresponding e. The result is achieved very fast (9 iteration). The 
best performance that can be achieved by new method is ^^ew — 0.8345 at e = 0.0802, which is 
about 36% improvement, over tli6 conventional metliod. The controller associcite with, 'y^eii? 

K = [ -1.2732 -0.8923 -1.8967 -1.8145 ] 



Searching for three independent e 's can be a tedious job, particularly if the number of 
objectives is rather large. Here, to show the effect of independent e 's, we fix one of the e's, 
starting with the value obtained above (i.e., when all were set equal to one another), and search 
for the other two e 's, which are assumed to be equal. The best result obtained is 7* = 0.8104 
for £2 = 0.0802 and ei = £3 = 0.1292. Searching for three independent e leads on a slightly 
better result (about 3%), which might be enhanced if a full - or more thorough - search is done, 
but the computation burden will be significant. 



3.3 Output Feedback Case 

In this section, we use the dilated Mis in designing a full order dynamic output feedback 
compensator for the same saturation problem. Therefore, the controller is 



(14) 



where Ac is of the same order as the system matrix A. To reduce clutter of the equations, we 
dropped Dc- Applying this controller, the closed-loop system in (3) can easily be obtained with 
X = [xl xlY. 

In output feedback case, additional complications arise since a variety of transformations 
and manipulations are needed to set the problem into a convex search, often requiring auxiliary 
variables instead of the compensator matrices {A^ Be and Cc). The approach is reasonably 
well known and can be found in a variety of references (e.g., [4] and among many). Here, 
most of the technical details are omitted but can be found in the references mentioned, though 
the outline is based on the approach used in [3] and p!3]. Here, for the Lyapunov matrix in 
Algorithm 1, we use the structure 



Y -Y 
-Y S-^ + Y 



> 



which, as discussed in [14], can be done without any loss of generahty. Therefore, Algorithm 1 
can be expressed through the following lemma. 



Lemma 5 ([2]) System ^ with disturbance w(t) satisfying w{t)'^w{t) < w^ax '^'^^ with a 
output feedback compensator is internally stable, never saturates and has a disturbance 
attenuation level jcon if there exist y > 0, X > 0, and general matrices L, F and E and a 
constant aeon > such that 
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where U = He{AX + B2Y) and A = He(A'^y + EC2). Then for S = X - Y"^ one represen- 
tation of the controller matrices is: 



Cr = FS 



-1 



Br 



-y-'E 



Ar 



{A - BcC2)xs-^ + B2CC - y-^LS-^ 



For dilated Mis, also we need to do some manipulations to expand them into appropriate 
convex or near convex forms. We start with dilated MI for L2 gain ([5]). As in the state feedback 
multi-objective solution, eventually we are going to use common G, the key slack variable 
introduced by dilation. As a result, for simplification, the index i of Gi {i = 1, 2, 3) is dropped. 
Let us call G~^ = H. Note that based on the structure of dilated Mis, G is invertible. By pre- 
and post-multiplying (UD by 

Diag [ T^H^ I I T^H^ ] 
and its transpose, respectively, we obtain 

/ $11 



T'^H'^XiHT + UeiT'^ H"^ AciT) - i$44 
T^{H^XiH + H- 2eiiH^A,i - iff^))T 
$44 = ne{T^H^T) 

Here, T is an auxiliary matrix used for the additional transformations that are needed in the 
output feedback synthesis problem. 

We partition H and G into the following forms with each submatrix having the dimension 

n X n 

H = ( I G = ( I (16) 

\ H21 H22 J ' \ G21 G22 J 



r'H^'Ba 
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(15) 



where 

$11 

$14 = 



As mentioned before, G is invertible, so the sub-matrices of this matrix are also invertible (by 
invoking a small perturbation if necessary [IS])- Now, to turn dilated Mis into near convex 
form, let's consider the the auxiliary transformation matrix T as 

rn _ ( Gil I \ 
\G2i )■ 

To simplify notations, we call Gu = R and Hn = Y. Finally, let us call 
T-HXiHT = M,= [f^l ^;^l);iMi = Mr) 

and 

F = C,G2i, E = hIb„ V = R^Y + Gl^H2i 

L = {Y^A + EC2)R + Y^B2F + H^^A^Gii 

These are the new variables introduced to turn this Mis into convex form (similar to G, -F, L 
used in the standard approach. Also to save space, for i = 1,2, we use 

11 = AR + B2F - \r, k = ya + ec2- \y^ 

n = + V) + A + , = -Mil + R^ - 2ein 

Ei = -Mi3 + Y - 2eiA, Ti = -Mi2 + V-2ei{A- h) 

Ti = -Ml + I-2ei{L-]^V^) 

Naturally, the same manipulation can be conducted on MI (jlOp and ()12p . albeit with a bit 
less clutter. Considering above matrices and definitions, the entries of the three Mis can be 
expanded in detail as in the following lemma which expresses Algorithm 2 for dynamic output 
feedback compensator. 



Lemma 6 System ^ with disturbance w{t) satisfying w'^{t)w{t) < w"^^^ and with a output 
feedback compensator ^14\ ) is internally stable, never saturates and has a disturbance attenuation 
level 'jnew if there exist square matrices R, S, V and symmetric matrices Y , Mi, M2 
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and general matrices L,E, F, and constant a. 
such that 

/ Mil + He(n) M12 



> 0, and small positive value < 1 (z = 1, 2, 3) 
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Then for invertible G21 and H21 deduced by matrix factorization of G21H21 = V — R^Y as 
mentioned in 116], the controller matrices are as follows: 



Cr — FGo-^ 



21 1 



H21E 



Ac = H^^iL - Y' AR + EC2R + YB2F)G2i 



-1 



As mentioned above, invertible G21 and i/21 can be deduced by matrix factorization of 
G21H21 = V — R^Y. As indicated in [16], this deduction is always possible and if necessary 
we can use perturbation. In our numerical examples, when possible, we obtain G21 and H21 by 
picking H21 = Y and therefore having G21 = Y~'^V'^ — R. 

As before, to avoid excessive computational cost in the search for ei, €2, and £3 in Lemma 
([6]), we can use same e for all Mis. 



3.4 Numerical example: Output Feedback 

Consider the same numerical example as for the state feedback, this time with k = 0.4, / = 0.04, 
uiim = 100 and with = [xi 3:2]. 

This system is expected to withstand disturbances with peak bound of Wmax = 5. Using 
the conventional approach, the resulting controller has Ac and Be matrices of order 10^ and 

Cc\con = [ -12.5082 -18.5711 -5.8239 -42.4745 ] 

The induced L2 gain of this system under this controller is 7*0^ = 1.7636. 

This problem is also solved with the new dilated matrix inequalities. We start with the 
same e in all three dilated matrix inequalities, and we use Golden-ratio method in search over e. 
This, after 10 iteration, leads to "^new — 1-5029 at e = 0.0231, which is about 15% improvement, 
over the conventional method. The corresponding controller has A^. and Be matrices of order 
10^ and 

Gc\new = [ -11-5857 -17.0741 -5.5596 -38.7530 ] 



Next, we try independently varying e's. To limit the cost of searching for three independent 
e's, we have fixed one of the e's in a certain value and search for the other two e's which 
are assumed to be equal. The best result obtained is 7* = 1.2746 for 62 = 0.0181 and ei = 
63 = 0.1618. Thus, letting e's vary independently has a significant effect at the results, an 
improvement of about 27% over the conventional method of solving the problem (7*0^ = 1.7636). 



4 CONCLUSIONS 

We presented new dilated matrix inequalities for Bounded Real MI, invariant set MI and con- 
straint MI. The structure of these dilated Mis, in which system matrices are separated from 
Lyapunov matrices, allows us to use different Lyapunov matrices for diflPerent objective in multi- 
objective problems or different parameter values in robust synthesis problems. The new ap- 
proach is guaranteed to achieve results which are better or equal to the ones obtained from 
the standard multi-objective setting. The synthesis results, for both state feedback and output 
feedback problems, are demonstrated through an example. 
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